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Abstract 

We study the probabihty distribution of the maximum Ms of a smooth sta- 
tionary Gaussian field defined on a fractal subset S of M"". Our main result is the 
equivalent of the asymptotic behavior of the tail of the distribution ¥(Ms > u) as 
u — )• -|-oo. The basic tool is Rice formula for the moments of the number of local 
maxima of a random field. 

1 Introduction 

Let {X{t),t G M} be a stationary centered Gaussian process with constant variance equal 
to 1 and paths and let M be its maximum on a compact interval, for example [0, 1]. 
It is well-known since the pioneering work of S.O. Rice [17] that 

P{M >«}<!- <I>(m) + (const) (f{u), 

where ip is the density of the standard Gaussian law, $ its cumulative distribution func- 
tion and (const) is a constant depending upon the law of the process. In what follows, 
(const) will denote a constant that may change in each occurrence. 

Under stronger hypotheses, one can prove the equivalence (see, for example, Piterbarg 
[H], Azais and Wschebor [6] Prop 4.1 and 4.2). 

P{M > u} ^ (const) (p(u) as m — > -|-oo. 
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We are using the following definitions and notation: For /, g real- valued functions of one 
real variable, we put {xq can be a real number, +00 or —00): 



f — g as x — )■ xo 



if lim. 



■x—^xo 



f{x)/g{x) = 1. 



Similarly, we put 



f ^ g as X — )■ Xq 



if there exists a positive constant c such that 1/c < f{x)/g{x) < c for all x 7^ Xq in some 
neighborhood of xq. 

Let now X = {X{t),t G R"}, n > 1 be a sufficiently smooth real-valued centered 
stationary Gaussian random field, with a law that satisfies certain non- degeneracy as- 
sumptions. We consider a compact set S, C M" having topological dimension d and 
possessing a geometry that satisfies a certain number of smoothness properties that we 
do not detail here. Let Ms be the maximum of the random field on 5*. It is known (See 
for example Adler Piterbarg [T6], Taylor, Takemura and Adler [18], Adler and Taylor 
^ Azais and Wschebor [5], |6j) that: 



Under certain conditions, it is possible to go beyond this equivalence and obtain a 
whole expansion for u — )■ +00. In fact, the approximation becomes very accurate, in the 
sense that for large values of u, the error term is exponentially smaller than the principal 
term. Other extensions are known, for example to non- stationary Gaussian fields ([2], [S], 



The goal of this paper is to extend formula ([T]) assuming that S* is a fractal set, in 
which case no results have been published yet, as far as the authors know. We will as- 
sume in the sequel that S" is a bounded subset of M" and denote M5 = sup{X(t) : t E S}. 
Assuming that the paths t -w X(t) are continuous, it is easy to see that Ms is actually a 
random variable, without any further hypothesis on S. 

We will prove that whenever 5* is Minkowski measurable in the strong sense (Falconer 
[TU] . definition below) formula ([T]) holds true replacing d by the Minkowski dimension 
and the constant factor can be explicitly computed from the law of the process and the 
content of the domain. This includes a large class of self-similar sets. In other cases, we 
are able to obtain the asymptotic behavior of f{Ms > m) as m — >■ -foo up to a constant 
factor or - more weakly - the equivalent of logP(M5 > u). 

Another interesting situation is when 5* which has an outer Minkowski content but 
may have a non-smooth boundary (see example 4 below). In this case, we obtain the 



P{M5 > m} ^ (const) u'^ Vl"") as m -> +00. 



(1) 



m)- 
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asymptotic behavior of F{Ms > u) for large m, which depends upon the Lebesgue mea- 
sure and the perimeter of the set S. 

We are only considering here equivalences of the form ([1]) not full expansions which 
appear to be essentially more complicated when the set S is fractal. This will be the 
object of future work. 

For each subset S of M" and each e > 0, we will denote S*^ the tube of radius e 
around 5, that is the set S*" = {t G M" : dist{t,S) < e}, where dist denotes Euclidean 
distance. This set is sometimes called "parallel set" or "e-neigbourhood". The set S 
is called "Minkowski measurable in the strong sense", with Minkowski dimension d and 
content C, if 

X^S') ~ Ce"-^ as e ^ 0, 
A„ denotes Lebesgue measure in M". 

Similarly, 5* is called "Minkowski measurable in the weak sense", with Minkowski 
dimension d if 

XniS') ~ e"-"^ as e ^ 0. 
More generally, the Minkowski dimension d is defined as: 

logA„(g-) 

d = n — lim , 

e-s>0 log e 

whenever the limit exists. 



2 Main theorem 

Let -B be a ball centered at the origin such that S G B. We set B = B^. 

For M G M consider the excursion set 

Eu = {teB: X{t) > u}. 

Let t be a local maximum of that is, 

X{t) > X{s) W s E B n U for some open ball U centered at t. 

If t G Eu, we denote by Kt the connected component of Eu that contains t. Note that 
such a t can lie on the interior B of the ball B or on its boundary dB. 

We make the following assumptions on the random field X : 

• It is Gaussian, centered and stationary. We denote its covariance by 
T{t) =E{X{s)X{s + t)). 
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• The paths are of class C^. 

• For aR s y^t & B the random vector: 

{X{s),X{t),X'{s),X'{t)) has a non degenerate distribution. 

• For alH G -B, AG 5"^^ the random vector: 

{X{0),X'{0),X"{0)X) has a non degenerate distribution. 

• With no loss of generality wc assume in addition that it is centered, and for every 
t one has Var(X(t)) = 1, Var(X'(t)) = /„, where /„ is the identity matrix. Con- 
cerning this last condition, in any case, one can make a linear transformation of the 
parameter space so that it is satisfied, changing the set S into its image under this 
transformation . 

If A is a unit vector in <S"-^ we denote = X^X"{s)X. B{t,r) {t eW, r > 0) 

denotes the open ball with center t and radius r in R". 

The following lemma is essential for the proof of the main theorem. 

Lemma 1 Let a be a real number < a < 1. Let be the number of local maxima 
above u in B. Under our assumptions, the following quantities are "negligible" meaning 
by that that they are o{u~^ip{u)) as u ^ +00; 

1. P(Ai) where Ai — {3 a local maximum in B with value >u-\-l}. 

2. E(M„(M„-1)). 

3. P(vl2) where A2 = {3 two or more local maxima in B with value > u}. 

4. P(A3) where 

As = {3 a local maximum t & B 

s — t 

such that u < X{t) <u + l, min{X^'(s) : s G B{t, p), A = ^ —^} < -X{t) - u"} 

where 

p = p[u) = u'^ with /3 > (1 - a)/2 
This means that the connected component associated to t is not too small. 

5. P(^4) where 

= {3 a local maximum t G B 

s — t 

such that u < X{t) < -u + 1, max{X'^{s) : s G B{t, p),X= r. rA > -X{t) + -u"} 

s t 

This means that the connected component associated to t is not too large. 
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Proof : 

1 is a direct consequence of the fact that E{Mu) < {const)u'^^^ip{u) 

2 is proved in |15] and in Delmas' PHD dissertation [8] under shghtly more restrictive 
conditions on the random field, but the proof goes through with the same arguments 
under our hypotheses. The result is also quoted in 

3 is a direct consequence of 2. 

The proof of 4 is very similar but simpler than that of 5. 

Proof of 5. Let M4 be the number of local maxima in the interior of B such that 
u < X{t) <u + l and 

max{X;'(s) : s G B{t, p), A = ^} ^ -^(^) + 

Then ^(^4) < E(M4). This expectation can be computed by means of a Rice formula 
(see [B], Chap. 6): 

E(M4) 

^ In ^ ^^"W^O Imax{X;'(s):seB(*,p),A=^}>-.+„c|X(t) =X,X'(t) =0) 

PX{t),X'(t){x,Q)\n{dt). 

1a denotes the indicator function of the set A, -< means negative definite, Pci-) is the 
probability density function of the random vector C, . 

Using the stationarity of the law of the random field and the form of the density, we 

get: 
E(M4) 

< {const) I <^(s)e(| det[X"(0)]| llmax{X;'(s);sGB{0,p),A=^}>-x+«c. |X(0) = X,X'(0) = o) rfx 

(2) 

We use the Schwarz inequality for the conditional expectation: 

1 2 

E(|det[X"(0)]| I^ax{X;'(s):seB(0,p),A=^}>-x+«'^|X(0) = x,X'(0) = 0^ 

< E(^| det[X"(0)]|2|X(0) = x,X'(0) = 0' 

X P(^max{X;'(s) : s G B(0,p),A = ^} > -x + m"|X(0) = x,X'(0) = o). 

(3) 

For the first factor in the right-hand side of ([3]), note that X(0) and X'(0) are independent 
as well as X"(0) and X'(0). Regression formulas imply that, under the condition. 
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X"(0) = -Xin + C, 

where the matrix C = X"{0) + X(0)/„ is independent of X(0),X'(0). Hence: 

E(^|det[X"(0)]p|X(0) =x,X'(0) = o) < (const) [x^" + l] . 

Let us turn to the second factor in the right-hand side of ([3]). Under the condition, we 
have: 

X';is) = Tlis)x + R,is), 

where R\{s) is a random field with parameter s and bounded variance and T denotes the 
covariance function of X which is of class C^. Using standard properties of the covariance 
and the continuity of its fourth derivative, we can write 

^'i{s) < -1 + {const) \\sf. 

Taking into account our choice for a and /3, it follows that for large enough u we get the 
following bound for the second factor: 

p(^max{i?A(s) : s G B{0,p)} > - {const)\\sf{u + l)j 

< P(^min{i?A(s) : s G B{0,p)} < -{const)u' 

< (const) exp [ — (const) u^"] . 



(4) 



The last inequality follows from the classical Landau-Shepp-Marcus-Fernique inequal- 
ity for the tails of the supremum of a.s. bounded Gaussian sequences (see for example 
[6], (2. 33)). Replacing in ([2]) we finish the proof of 5. ■ 



Theorem 1 Assume that the process X satisfies the above conditions. Then, 

1. If the parameter set S is Minkowski measurable in the strong sense, with dimension 
d and content C , one has, as u ^ +oo; 

>u)^ ^^7^r(l + in- d)/2) u'^-Mu). (5) 

2. If the parameter set S is Minkowski measurable in the weak sense, with dimension 
d, one has, as u ^ +oo: 

P(M5 > m) ~ u^-^^{u). (6) 

3. If the Minkowski dimension d of the parameter set S is well-defined then, as u ^ 
+oo: 

logP(M5 >u) = --u^ + {d-l + o(l)) logM. (7) 
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Proof : We denote 

4 

A=[jA, 

1=1 

According to the Lemma, as m — > +00 

P(A) = o{if{u)/u). 

If u is large enough, whenever the complement set A'-' occurs and the excursion set is 
non-empty, the connected component Kt associated to the local maximum t E B satisfies 

B{t,r) ClKt(ZB{t,f) (8) 

Where 

. = .;2;^.'"-". .= ,/?<^ 



X{t) + V U-U"' ' 

This means that for u sufficiently large, such that ^J2/{u — u°') < 1 the component 
associated to a local maximum belonging to dB does not touch S. This implies that 
maxima in dB do not have to be taken into account when considering the event 
{Ms > u}n A^. 

Using Lemma 1, on A'~^ the excursion set has at most one connected component 
associated to a local maximum in B. 

Let us prove (IHl). Notice first that if u is large enough, one has B{t, r), -B(t, f) C -B(t, p). 
We prove the first inclusion in ([8]), the proof of the second one is similar. 

Let s G B(t,r), s ^ t. Write the Taylor expansion 

Xis) = Xit)+^-\\s-trX';{9) 
where A = (s — t)/\\s — t\\ and 6 is some point in the linear segment [s,t]. So, 

X{s) > X{t) -^-\\s- tf (X(t) + > X{t) - §^^iX{t) + U-) = u. 
This proves that s E Kt. 
Now consider: 

F{Ms > u) = o{ip{u)/u) + P^{3t G B,X{.) has a local maximum at t, u < X{t) < u + 1} 

(9) 

Then, an instant refiection shows that we have the following upper and lower bounds 
for P(M5 > u): 

P(M5 > m) < o(^(m)/m) + P({3t G B,X{.) has a local maximum at t,M < X{t) < m + l,t G S'^}), 

P(M5 > u) > o(v?(u)/u) + P({3t G has a local maximum at t,M < < u + l,t G S-}). 

(10) 
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Denote 

Mu^u+i,r = & B : t E S"", X{.) has a local maximum ed. t,u < X{t) < u + 1}. 
The first inequality in ( !T0|) implies: 

P(M5 > n) < o((^(?z)/n) + E(A4,„+i,,). 
As for the second inequality in (ITOl) . we have: 

< E(M„,„+i,, ) - P({M„,,+i,, > 1) < ^E(il/4,„+i(M„,„+i - 1)) = o{^{u)/u) 
so that 

P(M5 > > o{ip{u)/u) + E(M„,„+i,J 

To finish the proof we study the behavior as m — t- +oo of K{Mu^u+i,r) and E(M„^„+i ,.)• 
Notice that r,f are random. Let us consider case 1., when the set S is Minkowski mea- 
surable. Cases 2 and 3 are entirely similar. 

Let us look at E,{Mu^u+i,f) ■ Our tool is again Rice formula (See [6J, Chap. 6) which in 
this case takes the form: 



^ dx J^E(^\det[X"{t)]\ l{x"ituo} l{tesn/x{t)=x,X'{t) = o)pxw,x'w(x, 0)A„(di) 



We have, for u large enough: 



/ = ^^ dxj^ l^tes^^yE(\det[X"{t)]\l{x"it)-<o}\x{t)=x,X'{t) = 



-(p{x)\n{dt) 



u+l 



A„(5^"*)E(|det[X"(0)]| l{x' 



(0)^0} 



(27r)-/2^ 

X(0) = x,X'(0) = o)(^(x)dx (12) 



where we have replaced the density by its value and used that the conditional expectation 
is independent of t, since the process is stationary, f* is defined as the value of f when 
replacing X{t) by x. Note that it is non-random. 

It is well known that the conditional expectation in the right-hand side of ( IT2|) is 
equivalent to x" as x — )■ +oo. (See for example p]). Use now that the set S is Minkowski 
measurable in the strong sense, with dimension d and content C. This easily implies that 
as M — )• -|-oo one has: 
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1 



u+1 



C 1 



x'^exp(-xV2)C 



u — 



-1 (n-d)/2 



dx 



u 



{n+d)/2 



u+l 



{x-uf''-^^'^ex-p{-x^/2)dx 



Now perform the change of variables x = u{l + y /u^) m the last integral. We get: 



C 



c 



d-1 



n — d 



v.(«)r(^ + i) 



where the last equivalence follows using Lebesgue's Theorem. 

For the lower bound of f{Ms > u), we use again Rice formula to express K{MuM+i,r) 
as an integral. The computation is entirely similar and we obtain the same equivalent as 
M — 7- +00. This ends the proof. | 



3 Examples 

1. Smooth manifolds 

Our first example is the well-known case in which S* is a compact smooth manifold 
without boundary having topological dimension d embedded in M". In this case, 5* 
is Minkowski measurable in the strong sense, with dimension d and content equal 
to 

a,iS)VoliB„_,), 

where ad{S) denotes the geometric measure of S and Vol{Bn-d) is the volume of 
the unit ball in R""*^. Replacing in ([5]) we get: 

This is a much weaker version of results in [18] or ^ 

2. Self-similar sets. 

Let fi-.W ^ R" {i = 1, .., k) be k similarities of R", that is, 

\\Mx)-f,{y)\\=Ci\\x-y\\ V x.y G R" 

and we assume that < q < 1 for i = 1, k, k > 2. Let S be the self-similar set 
associated to fi, f^, which means that 

k 

s = [jMS). 

1=1 
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We also assume that the /j's satisfy the "open set condition", that is there exists a 
non-empty bounded open set V such that 

k 

[jMv)cv 

1=1 

where the union is disjoint. 

Under these conditions, one can prove that the Hausdorff and Minkowski dimensions 
of S are equal, and their common value d is the unique solution of the Moran 
equation: 

k 

1=1 

So, if S is self-similar and satisfies the open set condition, one has the logarithmic 
expansion ([7]). 

Classical examples of self-similar sets satisfying these conditions are the homoge- 
neous Cantor sets, the Sierpinski gasket, the von Koch and modified von Koch 
curves. For a proof of the above statement and a number of other examples and 
extensions, see [TU] ). 

3. Special self-similar sets. 

For self-similar sets one can obtain Minkowski measurability in the strong sense, at 
the cost of adding a certain number of conditions. One can compute the dimen- 
sion and the content and then obtain the equivalence for the tail of the distribution 
of the maximum having the form ([5]) in the statement of the main theorem . We 
restrict ourselves to dimension n = 1. The formulas become more complicated for 
higher dimensions, but it is possible to extend the results, adding a certain number 
of conditions (see [9]). 

In dimension 1, the result is given in p!3] and [Tl|. We use the notation of the previ- 
ous example and assume the same hypotheses. Let the interval I be the convex hull 
of S. Then, for i ^ i', fi{I) and /«'(/) are disjoint, except possibly at the endpoints, 
so that / \ U^^ifi{I) is the union of k', k' < k — 1 disjoint intervals with lengths 

£i, 

The set S is said to be of lattice-type if the additive subgroup of the reals Yl^=i (log Cj)Z 
is discrete. Then, the result is that S is Minkowski-measurable in the strong sense, 
if and only if it is not of lattice-type, and in that case, the content is given by: 

ol-dV^fc' fid 

c = -^'"-^ ^ 

rf(l-rf)EllCflog(l/Q)' 

As an illustrative example, assume that k = 2. Then, S is not of lattice type if and 
only if log Ci / log C2 is irrational. 
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4 Sets with a finite perimeter 

Let S C M" be a closed set. Whenever the hmit 

lini = OM{S) 

e^O e 

exists, S* is said to have the "outer Minkowski content OM{S)." 

For a review of this subject and related ones, see for example [3j. Let us mention here 
some general examples of classes of sets having outer Minkowski content: 

1. S* is compact with Lipschitz boundary, in which case 

OM{S) = DG{S) = H'^-^dS) (13) 

where DG denotes De Giorgi's perimeter (see [7]) and H"" the usual Hausdorff 
measure of exponent a. 

2. Let S* C M" be a closed set and define: 

U{S) = {a G M" : 3 a unique x E S such that dist{a, S) = dist{a,x)} 

reaches) = inf sup{r > : B{a,r) C U{S)}. 

S is said to have "positive reach" if reach{S) > and in this case, S has an outer 
Minkowski content (see [12] )• Compact convex sets with non-empty interior have 
positive reach and in this case, ( IT^ is also verified, allowing to compute the outer 
Minkowski content as a perimeter of the Hausdorff measure of the boundary. 

3. The outer Minkowski content need not coincide with the Hausdorff measure of the 
boundary. More precisely, it is possible to construct compact sets S* C with 
non-empty interior, positive reach, hence having outer Minkowski content, but such 
that OM{S) > H'^~'^{dS) (see for example [3], example 1, where this is done for 
n = 2). 

We have the following new result: 
Theorem 2 Assume that S has outer Minkowski content. Then, as u +oo: 

FiMs >u) = US)^^ u-Mn) + ,(„,?)X!f2i)/2 ^-V(^) (1 + 0(1)) 

The proof follows the one of Theorem [1], with the only change that one needs a 
somewhat more precise result on the asymptotic behavior of the conditional expectation. 
In fact, one can use that (see for example [8]) 



E(|det[X"(0)]| I|x"(oHo} |X(0) = a;,X'(0) =0) =x" + 0(a; 
and this suffices to adapt the proof of the main theorem to this example. 



n-2\ 
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